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O Astrophysical black hole candidates are thought to be the Kerr black holes 

?^ predicted by General Relativity, but the actual nature of these objects has still 

bJ!) to be proven. The analysis of the electromagnetic radiation emitted by a geo- 
metrically thin and optically thick accretion disk around a black hole candidate 

^ can provide information about the geometry of the space-time around the com- 

0\ pact object and it can thus test the Kerr black hole hypothesis. In this paper, I 

^ present a code based on a ray-tracing approach and capable of computing some 

^ basic properties of thin accretion disks in space-times with deviations from the 

O Kerr background. The code can be used to fit current and future X-ray data 

C\l of stellar-mass black hole candidates and constrain possible deviations from the 
Kerr geometry in the spin parameter-deformation parameter plane. 

> 

^ Subject headings: accretion, accretion disks — black hole physics — general 

c5 relativity — X-rays: binaries 



1. Introduction 



The predictions of General Relativity have been confirmed by experiments in Earth's 



gravitational field (Williams et al. 2004 Everitt et al. 2011), by spacecraft missions in the 



Solar System (Bertotti et al. 2003), and by accurate radio observations of binary pulsars 



(Weisberg & Taylor 2005 Kramer et al. 2006) (for a general review, see e.g. Will (2005)). 



In all these environments, the gravitational field is weak, in the sense that one can write 
Qtt = 1 + 4' with 101 ^ 1. The validity of the theory in the regime of strong gravity, when the 
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approximation |0| ^ 1 breaks down, is instead still unexplored. The ideal laboratory to test 
strong gravitational fields is the space-time around astrophysical black hole (BH) candidates 



(Bambi 2011c). 



In 4-dimensional General Relativity, uncharged BHs are described by the Kerr solution 
and are completely specified by two parameters: the mass, M, and the spin angular mo- 



mentum, J (Carter 1971 Robinson 1975). A fundamental limit for a Kerr BH is the bound 
\a\ < M, where a = J/M is the spin parameteij^ For \a\ > M, there is no horizon, and the 
Kerr metric describes the gravitational field of a naked singularity, which is forbidden by 



the weak cosmic censorship conjecture (Penrose 1969). On the observational side, there are 



at least two classes of astrophysical BH candidates (for a review, see e.g. Narayan (2005)): 
stellar-mass compact objects in X-ray binary systems (M ^5 — 20 Mq), and super-massive 
bodies at the center of every normal galaxy (M ~ 10^ — 10^ Mq). The measurement of the 
mass of the two classes of objects is robust, because obtained by dynamical methods and 
without any assumption about the nature of these objects. However, this is basically the 
only solid information we have. We think they are the Kerr BHs predicted by General Rela- 
tivity because there is no alternative explanation in the framework of conventional physics: 



stellar-mass BH candidates are too heavy to be neutron star (Kalogera 1996), while at least 



some of the super-massive objects in galactic nuclei are too heavy, compact, and old to be 



clusters of non- luminous bodies (Maoz 1998). 



The Kerr-nature of astrophysical BH candidates can potentially be tested with already 
available X-ray data, by extending the two most popular techniques currently used by as- 
tronomers to estimate the spin parameter of these objects: the continuum-fitting method 



(Zhang et al. 1997; Li et al. 2005; McClintock et al. 2011) and the Ka-iron line analysis 



(Fabian et al. 1995; Reynolds & Nowak 2003 Miller 2007). With the continuum-fitting 



method, one studies the thermal spectrum of geometrically thin and optically thick accre- 
tion disks: under the assumption of Kerr background and with independent measurements 
of the mass of the object, its distance from us, and the inclination angle of the disk, it is 
possible to infer the spin parameter a and the mass accretion rate M. The technique can 
be applied only to stellar-mass BH candidates, as the disk's temperature goes like M~°'^^ 
and the spectrum turns out to be in the keV-range for objects in X-ray binary systems, 
and in the UV for the super-massive bodies at the centers of galaxies. Relaxing the Kerr 
BH hypothesis, one can investigate the geometry of the space-time around the BH candi- 



date (Bambi & Barausse 2011a). With the same spirit, one can analyze the broad Ka-iron 



lines observed in both stellar-mass and super-massive BH candidates to test their nature 



^Throughout the paper, I use units m which Gn = c = 1, unless stated otherwise. 



3 



(Johannsen & Psaltis 2012). Actually, one can usually get only a constraint on a certain 



combination of the spin parameter and of the deviations from the Kerr background. In order 
to really test the nature of the compact object, at least two independent measurements are 



necessary (Bambi 2012b). 



With this paper, I am going to present a code designed to compute the thermal emission 
of a geometrically thin and optically thick accretion disk around a compact object charac- 
terized by a mass, an arbitrary value of the spin parameter (that is, no subject to the bound 
\a\ < M, valid only for Kerr BHs (Bambi 2011a||b )), and by one (or more) "deformation 
parameter (s) " , measuring possible deviations from the Kerr geometry. The code is based 
on a ray-tracing method and takes all the relativistic effects into account, as well as some 
important astrophysical effects. It thus includes significant improvements with respect to 



the code used in Bambi & Barausse (2011a) and can analyze real X-ray data providing the 



most reliable test we can get with current knowledge. The analysis of the soft X-ray compo- 
nent of specific sources and the corresponding constraints on the spin parameter-deformation 



parameter plane will be presented in a forthcoming paper (Bambi et al. 2013). The code 



can also be used to compute other disk's properties, like its direct image and its light curve 
during an eclipse, which may be observed with future X-ray experiments and can provide 
additional details about the actual nature of a BH candidate. 



2. Thermal emission of thin accretion disks 



Geometrically thin and optically thick accretion disks are described by the Novikov- 
Thorne model ( Novikov fc Thor"iie||1973 Page fc Thorne|[l974 ), which is the relativistic gen- 



eralization of the Shakura-Sunyaev model (Shakura & Sunyaev 1973). Accretion is possible 



because viscous magnetic/turbulent stresses and radiation transport energy and angular mo- 
mentum outwards. The model assumes that the disk is on the equatorial plane and that the 
disk's gas moves on nearly geodesic circular orbits. The time-averaged energy flux emitted 



from the surface of the disk is (Page & Thorne 1974) 



M 



where F{r) is the dimensionless function 



F{r) 



(1) 



F r) 



/ {E -VLL){dpL)dp. 



(2) 



E, L, and fl are, respectively, the conserved specific energy, the conserved axial-component 
of the specific angular momentum, and the angular velocity for equatorial circular geodesies 
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(in a generic stationary and axisymmetric space-time, these quantities can be computed as 
described, for instance, in Appendix B of Bambi & Barausse (2011a)); G = —a^grrg^<i> is 



the determinant of the near equatorial plane metric, where = gt^/gcfxf, — gtt is the lapse 
function; is the inner radius of the accretion disk and is assumed to be the radius of the 
innermost stable circular orbit (ISCO). 

Since the disk is in thermal equilibrium, the emission is blackbody-like and we can 
define an effective temperature Teff(r) from the relation J-'(r) = aTjj, where a is the Stefan- 
Boltzmann constant. Actually, the disk's temperature near the inner edge of the disk can 
be high, up to ~ 10^ K for stellar-mass BH candidates, and non-thermal effects are non- 
negligible. That is usually taken into account by introducing the color factor (or hardening 
factor) /col. The color temperature is Tcoi(r) = fcoiTes and the local specific intensity of the 
radiation emitted by the disk is 

^e(^e) = .J . , (3) 

where z/g is the photon frequency, h is the Planck's constant, c is the speed of light, /cb is the 
Boltzmann constant, and T is a function of the angle between the wavevector of the photon 
emitted by the disk and the normal of the disk surface, say C,- The two most common options 
are T = 1 (isotropic emission) and T = ^ + ^ cos^ (limb-darkened emission). 



2.1. Calculation method 



The calculation of the thermal spectrum of a thin accretion disk has been extensively 



discussed in the literature; see e.g. Li et al. (2005) and references therein. The spectrum can 



be conveniently written in terms of the photon flux number density as measured by a distant 
observer, Ne^^^- A quite common approximation in the literature is to neglect the effect of 
light bending, which is good for small disk's inclination angles (face-on disks). In this case, 
Ne^.. is: 



Ne„ 



Eohs 

1 

Eohs 

A, 



5(^/e(z/o)c?fiobs 



-5'obs \ ' 

keV ) 



cos 2 



T\/—Gdrd(j) 



disk 



M2 



exp 



V kcV / 



(4) 



where /obs, -^obs, and u are, respectively, the specific intensity of the radiation, the photon en- 
ergy, and the photon frequency measured by the distant observer, cif^obs ~ cos i\/—Gdrd(j)/ 
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Geometry of the background 



Relativistic effects 

Self-irradiation 
Non-zero torque at 
Indirect images 
Color factor fco\ 
Radiation emission T 



Johannsen-Psaltis space-time with three parame- 
ters (mass M, spin parameter a, deformation pa- 
rameter €3). When £3 = 0, the geometry of the 
space-time reduces exactly to the Kerr solution. 
No restrictions on the values of a/M and 63 
All relativistic effects are included. Ray-tracing 
technique used 
Not included 
Not included 
Not included 

Constant. Set by the user 
Isotropic or limb-darkened 



Table 1: Basic features of the code to compute the observed spectrum of the disk. Indirect 
images are the images formed by null geodesies penetrating the equatorial plane inside rm. 
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1.2801 {9 = n/2) 
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1.2 


-1 


1.3614 {9 = n/2) 


2.5036 
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Table 2: Properties of the backgrounds shown in Figs. [IjjTj spin parameter a (first column), 
deformation parameter €3 (second column), radius of the event horizon th (third column), ra- 
dius of the ISCO risco (fourth column), and stability of the ISCO (MRS/MVS = Marginally 
Radially /Vertically Stable; fifth column). When 7^ 0, the radius of the event horizon in 
general depends on the angle 9 and the topology of the event horizon may be non-trivial (see 
Bambi fc Modesto] ( |2011| )). For the cases (a/M, €3) = (0.8, 2) and (0.9, 1), there is no event 



horizon on the equatorial plane 9 = 71-/2. 
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Fig. 1. — Observed thermal spectrum of a thin accretion disk in Kerr space-time. Left paneh 
a/M = 0, 0.5, 0.9, and 0.999 and viewing angle i = 30°. Right panel: a/M = 0.9 and viewing 
angle i = 0°, 40°, 70°, and 85°. The other parameters are M = 10 Mq, M = 10^^ g s'^ 
D = 10 kpc, /col = 1, and T = 1. These two plots should be compared with the two panels 
of Fig. 5 in Li et al. (2005). Flux density A^^^;^^^ in 7 keV~^ cm~^ s~^; photon energy i^obs in 
keV. 




Fig. 2. — Left panel: Kerr space-time with M = 1.0 ■ 10^^ g s ^ and 1.2 • 10^^ g s ^. The 
other parameters are M = 10 Mq, a/M = 0.999, i = 30°, D = 10 kpc, /coi = 1, and 
T = 1. The outer radius of the accretion disk is at rout = 10^ M. Right panel: Effect of 
light bending on the observed spectrum in Kerr space-time for a viewing angle i = 10° and 
80°. LB/NLB = light bending/no light bending. The other parameters are M = 10 Mq, 
a/M = 0.9, M = 1.0 ■ 10^^ g s-\ D = 10 kpc, /coi = 1, and T = 1. The outer radius of the 
accretion disk is at rout = 10^ M. Flux density Ne^^^^ in 7 keV~^ cm~^ s~^; photon energy 
Eobs in keV. 



-7- 



is the element of the sohd angle subtended by the image of the disk on the observer's sky, 
i is the viewing angle of the observer, and D is the distance of the source, g is the redshift 
factor 

where = hue, k" is the 4-momentum of the photon, u^y^^ = (—1,0,0,0) is the 4-velocity 
of the distant observer, and = 0, 0, I^m*) is the 4-velocity of the emitter. /e(z/e)/z/g = 
-^obs(i^obs)/i^^ follows from the Liouville's theorem. Ai and A2 are given by (for the sake of 
clarity, here I show explicitly Gn and c) 



Ai = 



2(keV)^ /GnM^ ^ 



0.07205 / M V /kpc 



2 

-1 ™-2 „-l 



kcV \ (G^mV''^ fAnaV'^ 



7 keV cm s 

/ kcV \ (G^MV'' UnaX 
\Mf^xJ \~J \W) 

1/2 

(6) 



Ao = 



/col \ M J 
Using the normalization condition Qi^^u'^u'^ = — 1, one finds 



0.1331 /10i«gs-i\'/' / M^'/' 



and therefore 



(7) 



^ = TTAf^ ' 

where A = /Ct^/Zct is a constant of the motion along the photon path. Neglecting the effect 
of light bending, A = r sin(/)sini, where r and arc the coordinates on the equatorial plane 
of the emitter. Doppler boosting, gravitational redshift, and frame dragging are entirely 
encoded in the redshift factor g. 

The effect of light bending can be included by using a ray-tracing approach. The initial 
conditions (io, ^^o; ^O; 0o) for the photon with Cartesian coordinates (X, Y) on the image plane 
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of the distant observer are given by (Johannsen & Psaltis 2010) 



to = 0, (9) 

ro = VX2 + Y^ + D\ (10) 

Y sini + D cos i 

On = arccos , , (11) 

00 = arctan . . — :. (12) 

D sm I — Y cos I 

As the initial 3-momentum kg must be perpendicular to the plane of the image of the 



observer, the initial conditions for the 4-momentum of the photon are (Johannsen & Psaltis 



2010) 



K = ^ =|ko|, (13) 

one — n y sini+Dcosi 

K = / ^x^+y^ ±Dl — (14) 

^X2 + (Dsinz-Fcosi)2' 
A + (iJ sm z — r COS «)"' 



kl = ^{k^o)' + rl{klY + rlsin'e,{kty. (16) 

In the numerical calculation, the observer is located at -D = 10^ M, which is far enough 
to assume that the background geometry is flat and therefore k^ can be inferred from the 
condition g^^k^k^ = with the metric tensor of a fiat space-time. The photon trajectory 
is numerically integrated backwards in time to the point of the photon emission on the 
accretion disk: in this way, we get the radial coordinate at which the photon was emitted 
and the angle ^ between the wavevector of the photon and the normal of the disk surface 
(necessary to compute T). Now g = g{X, Y), as in Eq. ^ everything depends on r^, except 
A that can be evaluated from the photon initial conditions (A = k^/kf = tq sin^o^o/^o)- 
The observer's sky is divided into a number of small elements and the ray-tracing procedure 
provides the observed flux density from each element; summing up all the elements, we get 
the total observed flux density of the disk. In the case of Kerr background, one can actually 
exploit the special properties of the Kerr solution and solve a simplified set of differential 
equations. That is not possible in a generic non-Kerr background and so the code solves 
the second-order photon geodesic equations of the space-time, by using the fourth-order 



Runge-Kutta-Nystrom method (Lund et al. 2009). The photon fiux number density is given 
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by 



Ne 



Eohs 

I Eohs 



-^obs('^)c?^obs 
2 



A, 



VkeV 



Eohs 

TdXdY 



g^Ie{l^c)dQohs 



M2 



exp 



( ^obs '\ 

V keV / 



(17) 



where X and Y are the coordinates of the position of the photon on the sky, as seen by the 
distant observer; that is, dQohs = dXdV/D'^. The basic features of the code are outhned in 
Tab.E 

The results of the code in the Kerr background are summarized in Figs. [T] and |2| while 
Tab. [2] shows some fundamental properties of these space-times. The two panels of Fig. [T] 



should be compared with the ones in Fig. 5 of Li et al. (2005). The agreement between 
the two sets of spectra is very good, except for the case a/M = 0.999 and inclination angle 
i = 30°, for which the difference is however small. The discrepancy is due to the fact that 



the present code does not include the effect of self-irradiation of the disk. As discussed in Li 



et al. (2005), the effect of self-irradiation of the disk, as well as a possible non-zero torque at 



the inner edge of the disk (also not included here, see Tab. [T]), can be ignored in the analysis 
of observational data, because they can be absorbed by adjusting the mass accretion rate 
and the spectral color factor of a zero torque model without returning radiation. This is 
indeed our case, as shown in the left panel of Fig. |2j where we can see the spectra of a Kerr 
BH with a/M = 0.999 and i = 30°, respectively for M = 1.0 ■ 10^^ and 1.2 ■ 10^^ g s-\ 
The spectrum with higher mass accretion rate and without self-irradiation of the disk is in 

(EOOSI) for M = 1.0 ■ 10^9 g s"^ and with self- 



Li et al. 



agreement with the one calculated in 
irradiationl^ The effect of light bending can be seen in the right panel in Fig. [2| which shows 
the observed spectra as computed respectively from Eq. ^ and Eq. (17). The difference in 
the observed spectra is significant for large inclination angles (edge-on disks) and negligible 
in the opposite case, when the angle between the observer and the symmetry axes of the 
system is small (face-on disks). In all these spectra, the inner edge of the disk is at the ISCO 
radius (Novikov-Thorne model), while the outer edge is assumed at rout = 10^ M. The latter 
is large enough that for the energy range shown in Figs. [TM] is equivalent to an infinite value. 



^With the standard approach in a Kerr background, one has two parameters, a and M, to be determined 



by fitting the soft X-ray component with the theoretical prediction of the thermal spectrum of a thin disk. Li 



et al. (2005 1 shows that disk's self-irradiation and small non-zero torque at the ISCO can be neglected without 



affecting the measurement of a: one just gets a slightly different M , which has not physical implications. 
As the physical mechanisms responsible to these effects are not peculiar of the Kerr geometry, the same 
conclusion should be true for metrics with deviations from Kerr. 
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The effect of Tout is indeed to introduce a different slope of the spectrum at low energies, see 



e.g. Fig. 4 in Bambi & Barausse (2011a). 



2.2. Non-Kerr space-times 



The Kerr-nature of astrophysical BH candidates can be tested by studying the thermal 
spectrum of the accretion disk without the assumption of the Kerr background. The basic 
idea is to consider a more general space-time, in which the compact object is characterized 
by a mass M, a spin parameter a, and a deformation parameter which measures possible 
deviation from the Kerr geometry. When the deformation parameter is set to zero, we 
have to recover exactly the Kerr solution. We can then compare the theoretical predictions 
with observational data: if the latter demand a zero deformation parameter, the Kerr BH 
hypothesis is verified; otherwise, we may conclude that the object is not a Kerr BH. 

Let us note that this approach can be used just to check if the geometry of the space- 
time is described by the Kerr solution, but it cannot really investigate the actual nature 
of the compact object or the exact deviations from the Kerr background. The analysis of 
the disk's spectrum can only probe the geometry of the space-time till the ISCO radius; we 
have no information about the geometry at smaller radii and about the surface/horizon of 
the BH candidate. Current (and near-future) data are not so good to map the space-time 
and a single deformation parameter is used to figure out if the gravitational force is stronger 
or weaker than the one around a Kerr BH with the same mass and spin. Actually, the 
typical situation is even worse and we can infer only one parameter: if we assume the Kerr 
background, we find the spin a, if we have also a deformation parameter, we constrain some 
combination of a and of the deformation parameter. However, this degeneracy can be solved 
with an additional measurement (Bambi 2012b|[cp ). 

When the computation of the properties of the thermal emission depends on the geom- 
etry of the background, the code calls the function metric. It is thus easy to change the 
space-time. For the time being, the code uses the Johannsen-Psaltis (JP) metric, which was 



explicitly proposed in Johannsen & Psaltis (2011) to test the geometry around BH candi 



dates. In Boyer-Lindquist coordinates, the metric is given by the line element (Johannsen 



& Psaltis 2011) 



2Mr 



[1 + h) dt^ 



AaMr sinH. S(l + /i) 

-(1 + h) dtd(j)+ \, . o ^ dr^ + 



+S de^ + 



sin^ 0\r^ + + 



2a^Mr sin^ 9 



+ 



a2(S + 2Mr)sin^^^' 
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where 



A 

h 



2. \ 2, 2 n 

r + a cos tl , 
- 2Mr + 



OO 

E 

fc=0 



Mr 

^2k H ?:r^2fc+l 



M2 

IT 



(19) 



This metric has an infinite number of deformation parameters and the Kerr solution is 
recovered when all the deformation parameters vanish. However, in order to reproduce the 
correct Newtonian limit, we have to impose eo = ei = 0, while €2 is strongly constrained by 
Solar System experiments. As done in other papers, it is usually sufficient to restrict the 
attention to the deformation parameter €3 and set to zero all the others. 

The properties of JP BHs with M, a, and €3 have been quite discussed in the recent 
literature ( |Bambi|[20lIdl |2012a||b1[cpl |Bambi k Modesto||201l| |Bambi et alT]|2012t [Johannsen 



& Psaltis 


2012 


Chen & Jing 


2012 


Liu et al. 


2012 


Krawczynski 


2012) 



these objects may have properties fundamentally different from the ones expected for Kerr 
BHs. For instance, when the gravitational force on the equatorial plane turns out to be 
stronger /weaker than the Kerr case (this corresponds to £3 < 0/e3 > 0), the maximum 
value of a/M can presumably be larger /smaller than 1, and increases/decreases as €3 de- 
creases/increases (Bambi 2011d). Another important feature to bear in mind is that the 



ISCO may not be determined by the orbital stability along the radial direction, as in Kerr, 
but it may be marginally vertically stable. The ISCO radius should thus be computed as 
described in Bambi & Barausse ( 2011a||b ). The event horizon of JP BHs has been discussed 
in Bambi & Modesto (2011) and the radius th can be found fror^ 



A + a^hsin^e = 



(20) 



For non-vanishing a and h, the radius th is a function of the angle 6. It is also possible to 
show that the topology of the event horizon of these BHs may be non-trivial and that BHs 



•^Let us note that there are a few definitions of horizon. Here we are talking about the event horizon, i.e. 
a boundary in the space-time beyond which events cannot affect an outside observer. In a stationary space- 
time, the event horizon is also an apparent horizon, which is a surface of zero expansion for a congruence of 
outgoing null geodesies orthogonal to the surface. This means that at the apparent horizon null geodesies 
must have dr/dt = 0, which implies g^^ = 0, see e.g. Poisson (2004 1. Eq. (20) follows from g^^ = in the 
case of the JP metric. The horizon relevant for the black hole thermodynamics is instead the Killing horizon, 
which is a null hyper-surface on which there is a null Killing vector field. For the metric in (18 1, the Killing 
horizon is defined by gttg<j><j> — gl^, — 0- When the Hawking's rigidity theorem can be applied (like in the Kerr 
space-time), the event horizon and the Killing horizon coincide ( Hawking|1972 1. However, in general that is 
not true. 
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with a topologically non-trivial event horizon may be created by astrophysical processes like 
the accretion from disk (Bambi & Modesto 2011[ ). The fundamental properties of the JP 
space-times shown in Figs. [3](7]are reported in Tab. |2} 

Examples of the observed spectrum of thin accretion disks in the JP background are 
shown in Fig. [sj for ajM = 0.8 and €3 = 2 (left panel) and for a/M = 1.2 and €3 = —1 
(right panel). Fig. [s] shows the effect of light bending in the JP background. As in the Kerr 
case, the effect is more and more important as the disk's inclination angle increases and the 
observer approaches the equatorial plane of the system. Fig. |4] shows instead the effect on 
the observed spectrum of the properties of the emission of the accretion disk. It compares 
the cases of isotropic (T = 1) and limb-darkened emission (T = ^ + | cos^, ^ being the angle 
between the wavevector of the photon emitted by the disk and the normal of the disk surface). 
The left panel of Fig. |4]can be compared with Fig. 9 of Li et al. (2005). As noticed in Li et 



al. ( 2005 ) , for edge-on disks the effect of limb-darkening cannot be absorbed by a redefinition 
of M and /coi; the spectrum has indeed a more pronounced hump before the exponential 
cut-off. The case of non-Kerr background does not introduce any new qualitatively different 
feature with respect to the Kerr case already discussed in the literature. 



3. Other observable features 



In the previous sections, I showed how the code can compute the thermal spectrum of 
a geometrically thin and optically thick accretion disk in a generic space-time, comparing 



the results for the Kerr case with the ones presented in Li et al. (2005). While even the 



continuum-fitting method is not immune to criticism, this approach to test the nature of 
astrophysical BH candidates is likely the best we can do with current observational data 



and theoretical knowledge (McClintock et al. 2011). Constraints on the spin parameter- 



deformation parameter plane obtained with the code by analyzing the data of specific sources 



in the high-soft state will be presented in a companion paper (Bambi et al. 2013). In this 



section, I am going to show some simple extensions of the code: the computation of the 
direct image of the disk and the prediction of the light curve during an eclipse. Current 
facilities cannot observe these features, but future X-ray experiments may do the job and 
provide additional information about the geometry of the space-time around BH candidates. 
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Fig. 3. — Effect of liglit bending on tlie observed spectrum in Johannsen-Psaltis space-time. 
Left panel: a/M = 0.8, = 2, and viewing angle i = 10° and 80°. Riglit panel: a/M = 1.2, 
es = —1, and viewing angle i = 10° and 80°. LB/NLB = liglit bending/no liglit bending. 
Tfie otlier parameters are M = 10 Mq, M = 10^^ g s~^, D = 10 kpc, /coi = 1, and T = 1. 
Tlie outer radius of tfie accretion disk is at rout = 10^ M. Flux density Ne^^^^ in 7 keV^^ 
cm~^ s~^; photon energy Eohs in keV. 




Fig. 4. — Effect of limb-darkened emission on the observed spectrum. Left panel: Kerr 
space-time with a/M = 0.999 and viewing angle i = 20° and 80°. Right panel: Johannsen- 
Psaltis space-time with a/M = 1.1, €3 = —0.5, and viewing angle i = 20° and 80°. The 
other parameters are M = 10 Mq, M = 10^^ g s'^, D = 10 kpc, and f^oi = 1-5. LDR/IR = 
limb-darkened/isotropic radiation. The outer radius of the accretion disk is at rout = 10^ M. 
The left panel can be compared with Fig. 9 of Li et al. (2005). Flux density Ne^^^^ in 7 keV^^ 
cm~^ s^^; photon energy Eohs in keV. 
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Fig. 5. — Direct image of the accretion disk. Observed blackbody temperature Tobs (left 
panels) and observed flux Jobs (right panels) in Kerr space-time with spin parameter a/M = 
(top panels) and 0.9 (bottom panels). The other parameters are M = 10 Mq, M = 10^* g s^^, 
i = 80°, and /coi = 1-6. The outer radius of the accretion disk is rout = 25 M. Tobs in keV; 
J-'obs in arbitrary units and logarithmic scale. 



- 15 - 




Fig. 6. — Direct image of the accretion disk. Observed blackbody temperature Tobs (left 
panels) and observed flux J-'obs (right panels) in Johannsen-Psaltis space-time with spin 
parameter a/M = 0.9 and deformation parameter = —1 (top panels) and 1 (bottom 
panels). The other parameters are M = 10 Mq, M = 10^^ g s~^, i = 80°, and fco\ = 1.6. 
The outer radius of the accretion disk is rout = 25 M. Tobs in keV; J-'obs in arbitrary units 
and logarithmic scale. 
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3.1. Direct image of thin accretion disks 



In the case of a geometricaUy thick and opticaUy thin accretion disk around a BH 
candidate, it should be possible to observe the so-called BH shadow; that is, a dark area 
over a brighter background. As the exact shape of the shadow is determined exclusively by 
the geometry of the space-time, its detection can be used to investigate the nature of the 



BH candidate ( 


Bambi & Freese 


2009 


Bambi & Yoshida 


2010 


Johannsen & Psaltis 


2010 


Bambi et al.||2012 


). In the case of geometrically thin and optically thick accretion disks, the 



image turns out to be quite different. For Schwarzschild and Kerr BHs, it has been already 



discussed in the literature (Luminet 1979; Fukue & Yokoyama 1988 Fukue 2003 Takahashi 



2004). The generalization to non-Kerr space-times is straightforward (see also Krawczynski 



(2012)). The temperature and the photon flux as detected by a distant observer are given 



by 



^obs 
•F abs 



(21) 
(22) 



where g is the redshift factor given by Eq. (|5|. Eq. (21) follows from the deflnition of g 
and from Eq. Q, while Eq. ^ is a consequence of the Liouville's theorem /e(t'e)/t'e — 
lohsiyohs)/^^- Fig. [5] shows the direct image of Tobs and J-'obs for a Schwarzcshild BH (top 
panels) and for a Kerr BH with a/M = 0.9 (bottom panels), in the case of an observer with 
inclination angle of 80°. The asymmetry with respect to the X axes is due to the effect of 
light bending, while the one with respect to the Y axes comes from the Doppler boosting. 
Fig. [6] shows the same images in the case the accretion disk is in the JP space-time with spin 
a/M = 0.9 and deformation parameter 63 = —1 (top panels) and 1 (bottom panels). As we 
can see, the shape of the central dark area and the intensity map of the image depend clearly 
on the geometry of the space-time around the compact object. Such images are deflnitively 
impossible to observe with current X-ray facilities, but they may be seen with future X-ray 
interferometry techniques. 



3.2. Light curve during an eclipse 



The observation of the light curve when the accretion disk around a BH candidate is 
eclipsed by the stellar companion can also provide information about the geometry of the 
space-time around the compact object (Fukue & Yokoyama 1988). While current observa- 
tional facilities would already have the correct time resolution to observe this feature, so 
far we know only one binary system with a BH candidate showing an eclipse in the X-ray 
spectrum. This source is M33 X-7, which is unfortunately quite far and dim, and present 
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detectors cannot measure its light curve with an accuracy to observe any relativistic effect 



(Pietsch et aL 2004). 



Fig. [7] shows the hght curves of the ingress (left panel) and egress (right panel) during 
an eclipse for the four cases shown in Figs. |5] and [6j The inclination angle of the observer is 
still i = 80°. In these simulations, I assume that the angular momentum of the companion 
star is parallel to the one of the accretion disk, as it should be more likely expected from 
considerations on the evolution of the system. That means that the companion star occults 
first the approaching (blue-shifted) part of the disk and then the receding (red-shifted) one. 
The companion star is simply modeled as an obstacle with vertical edge and its atmosphere 
is completely neglected, so the disk's luminosity at the time t is 

/rjXrIY 
h,{X-v,t)hU^)^^, (23) 

where i = ingress or regress, Vc is the velocity of the stellar companion, and 

,^ , r ifX-t;et<0, 
Wss(X-..t) = | ^ ifX-..t>0. 

^--^^-"=') = \ ifX-..t>0. 

As the shape of the light curve depends on the background metric, the possible detection 
of this feature can be used to constrain possible deviations from the Kerr nature of a BH 
candidate. However, in order to extract information on the space-time geometry from real 



data, a reliable atmospheric model would be necessary (Takahashi &; Watarai 2007). 



If we could get accurate images of the disk, like the ones shown in Figs |5] and [6| we 
could determine at the same time a and €3. However, that is surely out of reach in the near 
future. In the case of the light curve during an X-ray eclipse, basically we measure the slope 
of the curve and the asymmetry between ingress and egress. Assuming a Kerr background, 
we could immediately estimate a (assuming M is known independently). If we are going to 
test the Kerr-nature of the BH candidate, we find a degeneracy between a and €3. However, 
the light curve is essentially sensitive only to the effect of Doppler boosting; that is, the 
asymmetry of the apparent image of the disk with respect to the X = axes. It is not very 
sensitive to the effect of light bending (responsible to the asymmetry with respect to the 
y = axes) and to the gravitational redshift (which produces corrections symmetric with 
respect to the X = axes). The combination of the fit of the thermal spectrum and of the 
light curve during an eclipse for the same BH candidate should thus break the degeneracy 
between a and €3 and allows for the identification of a limited allowed region in the spin 
parameter-deformation parameter plane. 
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4. Conclusions 



General Relativity has been tested and verified in Earth's gravitational field, in the 
Solar System, and by studying the motion of binary pulsars. Thanks to recent theoretical 
progresses and high-quality data from present and near-future observational facilities, it is 
now possible to start testing the theory in the strong field regime. One of the most intriguing 
predictions of General Relativity is that the final product of the gravitational collapse is a 
Kerr black hole: a very simple object completely characterized by only two quantities (the 
mass M and the spin angular momentum J) in a very specific way. The study of the thermal 
spectrum of geometrically thin accretion disks around stellar-mass black hole candidates can 
provide information about the geometry of the space-time and can thus be used to check the 
Kerr black hole paradigm. 

In this paper, I have presented a code based on a ray-tracing approach and designed 
to compute the thermal emission of thin accretion disks around generic compact objects. 
In particular, the code can compute the thermal spectrum of a thin disk around a compact 
object with mass M, spin parameter a, and a deformation parameter which measures possible 
deviations from the Kerr background. By comparing these theoretical predictions with X-ray 
data of stellar-mass BH candidates in the high-soft state, one can constrain the nature of 
the compact object in the spin parameter-deformation parameter plane. Constraints from 



specific sources will be presented in a forthcoming paper (Bambi et al. 2013). The code can 
also be used to compute other observational features of a thin disk, like its direct image 
(Figs, [s] and |6]) or its light curve during an eclipse (Fig. [?]), which may be observed with 
future X-ray facilities. 



I would like to thank Luca Maccione for fixing a bug in the code. This work was 
supported by the Humboldt Foundation, Fudan University, and the Thousand Young Talents 
Program. 
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Fig. 7. — Light curves of ingress (left panel) and egress (right panel) during an eclipse for 
the cases of Figs. |5]and|6} Time t in units of M/vc-, where M is the mass of the black hole 
candidate and Vc is the velocity of the stellar companion. See the text for details. 



